Nuclear response for the Skyrme effective interaction with zero-range tensor terms. 

III. Neutron matter and neutrino propagation 
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The formalism of the linear response for the Skyrme energy density functional including tensor 
terms derived in articles [l|, |3| for nuclear matter is applied here to the case of pure neutron matter. 
As in article 0] we present analytical results for the response function in all channels, the Landau 
parameters and the odd-power sum rules. Special emphasis is given to the inverse energy weighted 
sum rule because it can be used to detect non physical instabilities. Typical examples are discussed 
and numerical results shown. Moreover, as a direct application, neutrino propagation in neutron 
matter is investigated through its neutrino mean free path at zero temperature. This quantity turns 
out to be very sensitive to the tensor terms of the Skyrme energy density functional. 

PACS numbers: 21.30.Fe 21.60.Jz 21.65.Cd 26.60.Kp 



I. INTRODUCTION 



In a recent series of articles P, , hereafter denoted 
respectively article I and II, the contribution of the zero- 
range tensor terms in the Skyrme effective interaction 
has been analyzed in the context of the Linear Response 
theory. The first result from these articles is that the ten- 
sor terms have very sizable effects on the response func- 
tions. Another important result is that the inverse energy 
weighted sum rule can be used as a tool of diagnosis for 
instabilities. These two articles were devoted to Symmet- 
ric Nuclear Matter (SNM) only. Since the construction 
of an Energy Density Functional (EDF) reliable for both 
symmetric matter and neutron matter is of fundamental 
importance 0, 01 , we present here the response functions 
and some associated sum rules for Pure Neutron Mat- 
ter (PNM) with the same approach. The interest of the 
present study is related to spin susceptibilities and fer- 
romagnetic finite size instabilities in neutron matter 0- 
[20j . Moreover, we use these results to study the impact 
of the tensor terms on the determination of the neutrino 
mean free path in PNM. This is a quantity of intrinsic 
importance since the cooling of a neutron star core in its 
first moments is governed by neutrino emission and there- 
fore by their mean free path trough dense matter. Some 
previous studies using non-relativistic approaches f2l| - [26j 
have revealed some very interesting features of the mean 
free path properties but they all neglected the possible 
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tensor contribution. Since the neutrino mean free path 
is directly related to the response functions which are 
themselves affected by the tensor, it is worthwhile to de- 
termine precisely the induced modifications. 

The article is organized as follows: in the first part de- 
voted to the linear response theory approach, we present 
explicit expressions for the spin response functions, the 
Landau parameters and the sum rules Mi, M3, M-i. 
Since the technical approach follows closely those of the 
previous articles, this part mainly contains figures and 
discussion, formula being written in appendices. The sec- 
ond part deals with the problem of neutrino mean free 
path. We first give explicit expression of this quantity in 
presence of tensor interaction then we show the influence 
of the parameterizations of the Skyrme functional. 



II. LINEAR RESPONSE APPROACH TO 
NEUTRON MATTER 



A. Response function 



Following article II, the starting point for the deter- 
mination of the response functions is the Skyrme en- 
ergy functional. Since in neutron matter the isospin is 
no longer a relevant quantum number and isovector and 
isoscalar densities are equal, it is convenient to define 
new coupling constants ~ Cq + Cf , where x ~ p, t, 
Ap, ... in such a way that the Energy Density Functional 
can be written as 
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with 
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The expressions of the coupling constants as functions 
of the parameters of the Skyrme interaction can be found 
in article I. From this expression, it is straightforward to 
obtain the residual interaction (see table in Appendix [5| 
by taking the second derivative of the EDF with respect 
to the density. 

The Random Phase Approximation (RPA) response 
function in each channel (a) = (S,M) = (spin, projection 
of the spin) (see Appendix [B|) is then obtained by solv- 
ing the Bethe-Salpeter equations for the correlated Green 
functions G^p^- Finally, from these response functions, 
we can easily obtain the Landau parameters (see Ap- 
pendix [C| . The quantities of interest are not directly the 
RPA propagators themselves, but merely the response 
functions S'*^"-'(q, w), also called the dynamical structure 
functions by some authors, which are defined at zero tem- 
perature by 



(3) 



From now on we choose the direction of q along the 
z axis, as done in article I and II. We show this func- 
tion in Fig. [1] for two different values of the transferred 
momentum {q ~ 0.05 fm"-'^ and q ~ 0.5 fm~^ and two 
different densities (p = 0.08 fm~'^ and p = 0.16 fm~'^) 
as a function of the energy w. As in article I and II we 
use a system of natural units so that h = c = 1. We 
consider one interaction without tensor, i.e. SLy5, and 
one with tensor, T16 (see article of Lesinski et al. [231 fo'^ 
the definition of the TIJ parameterizations) . Among the 
several TIJ possibilities the choice of T16 is motivated by 
the study of the neutrino mean free path (see later). In 
order to illustrate the effect of the interaction and of the 
RPA correlations we plot in each panel the corresponding 
Fermi gas (FG) and Hartree-Fock results (i.e. uncorre- 
lated response functions) . A first effect of the interaction 
clearly appears at the Hartree-Fock level where the mean 
field is responsible for the dressing of the bare neutron 
mass giving a density-dependent effective one. The differ- 
ence between the Fermi gas and Hartree-Fock structure 
functions increases with the density. With the RPA cor- 
relations, the difference between the S = Q and 5=1 
channel of the p-h interaction are reflected in the corre- 
sponding response functions. Let us focus on the 8=1 
channel, particularly important for the neutrino mean 
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FIG. 1; (Color on line) For the three PNM channels, we show 
the response functions S^°'^(q,oj) for the interaction SLy5 and 
T16. The vertical lines represent the position of the eventual 
zero-sound mode. For each force we plot as a reference the 
Fermi gas response function (dashed line) and the uncorre- 
lated response function (dashed-dotted line). 



free path. The {S = 1, M = 0) and {S = 1, M = 1) 
structure functions practically coincide, as it is expected 
for the SLy5 force, while they clearly differ in the pres- 
ence of tensor interaction. Concerning the 5 = 1 for 
low q (as illustrated in the figure for q = 0.05 fm~^) 
a spin zero-sound mode appears. It stands out above 
the p-h continuum. The existence of this spin collec- 
tive mode, called magnon, makes harder the excitation 
of the system, hence, correspondingly the p-h continuum 
response is depleted. This anti-ferromagnetic behaviour 
disappears when q increases. For high q another kind 
of divergence may appear. As illustrated in Fig. ^ the 
enhancement of the response function may become dra- 
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FIG. 4: (Color on line) Critical densities (in fm~^) as func- 
tions of the transferred momentum q (in fm~^) for the T16 
tensor parameterization and for SNM (panel (a)) and for 
PNM (panel (b)). The horizontal dashed-dotted line is placed 
at p = 0.16 fm~^ just to guide the eye. 



FIG. 2: (Colors online) The response functions ^'"'(a^) 
calculated for the Skyrme tensor parameterization T16 [27| . 
for the channel S = 1, M = only. The transfer momen- 
tum is g = 1.051 fm~^ and the density of the system is 
p = 0.16 fm-^ 
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FIG. 3: (Color on line) lEWSR (in MeV"^) as function of 
the transferred momentum q (in fm^^) for the T16 tensor 
parameterization. Full black line shows the result of the inte- 
gral Eq. Q while the dashed red line shows the result of the 
analytical expression Eqs. (|12p . (|13[1 and (|14p . The results are 
obtained at p = 0.16 fm~'^. The arrow indicates the position 
of the instability. 



matic and show a pole at a; = 0. In this case the homoge- 
neous Hartree-Fock ground state becomes unstable. For 
lower values of q the same kind of instability appears at 
higher density, called critical density pc, as illustrated in 
the next section. 
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FIG. 5: (Color on line) For the two 5=1 channels in PNM, 
critical densities (in fm~^) are plotted as functions of the 
transferred momentum q (in fm~^) for the TIJ family of 
Skyrme forces. The horizontal dashed-dotted line is placed 
at p = 0.16 fm"'' just to guide the eye. 



ticular the Mi Energy Weighted Sum Rule (EWSR), the 
M3 Cubic Energy Weighted Sum Rule (CEWSR) and 
the M_i Inverse Energy Weighted Sum Rule (lEWSR) 
defined as 



(4) 



B. Sum rules and moments of the strength function 

Following the notations of article II, we can calculate 
the most relevant odd-power sum rules for PNM, in par- 



As stated previously, all the expressions given below 
are derived for the general Skyrme EDF given by Eq.([3]) 
in which all the coupling constants could be considered 
as independent ones from the others. We refer to article 
II for a detailed discussion on their derivation. 
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The EWSR in each channel reads 
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Taking into account the M^2°'' value as well as the neutron 
effective mass m*, defined as 
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FIG. 6: (Color on line) For the three PNM channels, critical 
densities (in fm~^) are plotted as functions of the transferred 
momentum q (in fm"'^) for some usual Skyrme EDF: SkP 
SkM* [H], SGII [Ml, SLy4 [Hill.BSkS H and SkO 
The horizontal dashed-dotted line is placed at p = 0.16 fm 
just to guide the eye. 



— = 1 + 2m(CJ + Cl)p, 
m* 



(8) 



the Mf'^^ EWSR reduces to the free value /2m, as it 
should. As in article II for the case of a Skyrme force, 
these Ml moments can be obtained from the double com- 
mutator method [H, [s^l • 

For the CEWSR we have 
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And finally for the lEWSR we have 
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where the W^^^ 2 coefficients arc given in Appendix [B] 
As in article II we define the function /(fc) as 



f{k) = - 



l(l-fc^)lnf^ 
2fc ^ ' \k-l 



(15) 



while p and k are now defined for PNM with Fermi mo- 
mentum kp, hence 
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As already illustrated in the previous section, the main 
effect of the tensor terms is in the 5 = 1 channel where 
one can even observe a divergence at zero energy, but 
finite transferred momentum. As explained in detail in 
the article II, the lEWSR, M_i, can be used to detect 
these poles. As an example, on Fig. [3l we plotted the 
lEWSR for T16 obtained from the analytic expansion of 
the response function (see Eqs. ([T^ - P^ ') and from the 
direct numeric integration. We observe on this figure 



that in the channels S ^ 1, M = and 5 = 1, 71/ = 1 
the lEWSR is violated. This indicates the presence of a 
pole in the response function as shown for example for 
5^^''^^(g,w) on figurc[2]for p = 0.16 fm^"^ and transferred 
momentum q = 1.051 fm^^. 

This connection between the pole (when it does exist) 
observed in the response function and the pole observed 
in the M_i sum rule has been discussed in article II and 
we refer to it for a more detailed discussion on this point. 
It is thus possible to determine in a systematic way the 
critical density at which a pole occurs for a given mo- 
mentum q from Eqs. (fT2)) - (fT4l) . On Fig.[4]we display such 
critical densities, pc with respect to the transferred mo- 
mentum for the T16 interaction. On the left panel we 
first show the position of the poles of the response func- 
tion for SNM for each (S,M,I) channel. On the right 
panel we then show the position of the poles for PNM 
for each (S, M) channel. 

Even if we exclude the case of spinodal instability 
which is not present in PNM, one can see that the pres- 
ence and the location of the poles depends strongly on 
the system under analysis: for a given interaction, the 
critical densities are very different for PNM and SNM. 
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Similarly we show in Fig. [5] the critical densities for the 
tensor parameterizations that we use in the following sec- 
tion to study neutrino mean free path. 

Following article II, we show for completeness in Fig. [5] 
the critical densities for the Skyrme EDF previously an- 
alyzed in article II for SNM, but in this case for PNM. 
We observe that SkP behaves very differently from the 
other Skyrme forces presented in this article. It presents 
a first instability in the S — channel at p « 0.16 fm~^ 
and a second one at higher density p w 0.45 fm~'^ due 
to the presence of a pole in the effective mass defined in 
Eq. ®. 



III. NEUTRINO MEAN FREE PATH 



When tensor forces are considered the spin response is 
splitted into two components, i.e. the spin longitudinal 
response 



np 

and the spin transverse response 



np 



(19) 



(20) 



These responses can be considerably different one of the 
other and compared to the case without tensor interac- 
tion. This has important consequences on the neutrino 
cross sections and the neutrino mean free path. 



The aim of this section is to investigate the effect of the 
choice of the parameters of the tensor terms on the neu- 
trino mean free path in neutron matter. In this article we 
restrict ourselves to the academic case of neutron matter 
at zero temperature (the generalisation at finite temper- 
ature T is in progress). This, of course, implies some 
restriction on the application of our approach to neutron 
stars studies. At first stages of the cooling process of the 
neutron stars, modelized as asymmetric nuclear matter, 
the high temperatures involved allow charged current re- 
actions; then, when the temperature decreases, neutral 
currents dominate. Since we consider the case T = and 
the pure neutron matter, we will take into account only 
neutral currents for the determination of the neutrino 
mean free path. This quantity is defined as 

X={<7p)-\ (16) 

where a is the total cross section for the neutral current 
reaction v -\- n — > v' + n' . In the absence of tensor (and 
spin-orbit) interaction this total cross section is obtained 
by integrating the double differential cross-section per 
neutron 
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FIG. 7: (Color online) Axial spin transverse (AT), axial spin 
longitudinal (AL) and vector (V) contributions to the cross 
section per neutron for the reaction v + n — >■ u' + n' in 
the neutron matter. Three different densities as well three 
different interactions are considered. The Fermi gas result is 
also shown. 
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where Gp is the weak coupling constant, qa = 1.255 [s^ 
is the axial charge of the nucleon, Ei, and E'^, are the in- 
coming and outgoing neutrino energies, 9 the scattering 
angle between the incoming and outgoing neutrino mo- 
menta and uj = E,y — El and q = k — k' the energy and 
the momentum transfer in the reaction. The response 
functions R^{q,uj) and R^{q,uj) describe the response of 
the system to density fluctuations {S = 0) and spin fluc- 
tuations {S = 1) related to the coupling of neutrino to 
vector and axial currents of the neutron. These response 
functions are defined as 
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In the presence of tensor interaction the double differ- 
ential cross-section per neutron for neutral current reac- 
tion is given by 
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This expression reduces to Eq. (|17p when Rf^ = R^ = R^ 
as one can easily observe remembering that for neutral 
current 



q' = (k-k')^ 
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In Eq. ([2T|) . as well as in Eq. (fT7|) . we neglect cor- 
rections of order — from weak magnetism and other ef- 
fects [39| like the finite size of the nucleon or nuclcon 
excitations 5oj . A generalization of Eq. ([?!]) taking into 
account of all these effects can be found for example in 
Ref. . As already stressed in Rcf . (and illustrated 
in Ref. [4l| for charged current reaction) the cross sec- 
tion is dominated by the spin transverse response R^. In 
Fig. [7] we present the relative axial spin transverse, axial 
spin longitudinal and vector contributions to the neutral 
current cross section. We consider four different cases. 
The first is the Fermi gas. In this case = R^ = i?"* 
so the difference between the three contributions is only 
due to the kinematical factors and the coupling con- 
stants multiplying the response functions. Second we 
consider the case when the response functions are cal- 
culated with the SLy5 force. In this case the coupling 
constants = C^'^ = 0, purely related with the ten- 
sor part of the interaction, do not contribute. A possible 
difference between i?^ and Ri^ is due to the spin orbit 
contribution and according to the Egs lBSI and IB 71 given 
in Appendix [Bl the factor reduces these differences 
at low momentum transfer. Finally we consider the T13 
and T16 parameterizations of tensor contributions. The 
results for three different densities, p ~ 0.08, 0.16 and 
0.24 fm~'^ are also shown on Fig. [71 As it clearly appears 
all the cross sections are dominated by the spin trans- 
verse response but this contribution may vary between 

60 % and ^ 90 % depending on the interactions and 
densities considered. It reflects the possible quenching, 
enhancement or divergence of the nuclear response func- 
tions. This behavior was already illustrated in the section 
III Al in connection with the Fig. [1] To complete our dis- 
cussion we plot in Fig. [8] the spin transverse response 
for q = 0.05 fm"! and q = 0.5 fm^\ for p = 0.16 fm"^ 
and p = 0.24 fm'^ and for the tensor forces T13, T43 
and T63. A collective spin zero sound characterizes the 
S = I, M ^ ±1 response for T13 at g = 0.05 fm"! for 
p = 0.16 fm"'^ and p ~ 0.24 fm"''. A similar behaviour, 
with the corresponding quenching of the p-h continuum 
characterizes the 5 = 1, M = ±1 responses for SLy5 
and and T16 as shown on Fig. [T] for p = 0.08 fm~'^ and 
p = 0.16 fm~^. The T63 force on the other hand is char- 
acterized by an enhancement of the response at low for 
p = 0.16 fm~'^. At q = 0.5 fm~^ this enhancement seems 
to be critical. At p = 0.24 fm""^ the enhancement of the 
T16 response no longer holds. In this case this response 
is suppressed with respect to the corresponding Hartree- 
Fock response. An enhancement with respect to the HE 
and EG case for small w characterizes for this density the 
T43 force. 

These different behaviors affect obviously the neutrino 
mean free path. We have calculated it for T11-T16 and 
T13-T63 chains of parameterizations of Skyrme tensor 
interactions for a neutrino energy 5 MeV < < AO 
MeV and for three values of densities, i.e. p = 0.08, 0.16 
and 0.24 fm""^. Note that when a spin zero-sound col- 
lective mode appears one must in principle include it in 
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FIG. 8: (Color on line) The response functions for the channel 
5 = 1, M = 1 for three different forces (T13, T43 and T63) is 
shown with solid lines. On the same figure we represent the 
Fermi gas (dashed line) and the uncorrelated response func- 
tions (dashed dotted), i.e. when the residual interaction is put 
to zero, for the interaction T13. The uncorrelated response 
functions are not equivalents among T13, T43 and T63 due 
to the small differences in the effective mass. 



the calculation of the neutrino mean free path. For this 
mode the response function reduces to a delta distribu- 
tion. Nevertheless, as already observed in [2l|, the collec- 
tive mode itself gives little scattering, its contribution is 
negligible when calculated in the Landau approximation. 
It was also observed in [2^ that for all the Skyrme forces 
there considered, this magnon rapidly disappears with 
the temperature because of a strong Landau damping. 
Hence we do not explicitly compute the spin zero sound 
contribution. Its effect is on the other hand present as a 
suppression of the corresponding p-h continuum and has 
a consequence on the cross section. 

The results for the neutrino mean free path arc re- 
ported in Fig. [9] for the T11-T16 chain and In Fig. [10] 
for the T13-T63 chain. In each figure we include the 
Fermi gas result and the calculation with the SLy5 force 
which does not have tensor terms. From Figs. [51 and [TUj 
clearly emerges the crucial dependence of the neutrino 
mean free path on the parameterizations of the Skyrme 
tensor terms. For p = 0.08 fm~'^ the different parameter- 
izations give quite similar results and higher with respect 
to the non interacting case. Already at p = 0.16 fm~'^ 
the spread becomes important. For some cases A is higher 
than the Fermi gas one, for other is lower. In many cases 
it is lower than the corresponding SLy5 result. Also the 
behavior of A with the density is not so trivial, as al- 
ready appears in the three panels of Figs. 151 and [TUl This 
is clearly shown on FigITT] where the ratio A/ApQ is plot- 
ted. For example for T13 this quantity stays quite con- 
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FIG. 9: (Color online) Neutrino mean free path for scattering 
reaction u + n — > v' + n' in neutron matter. The dotted 
lines correspond to the non interacting Fermi gas case. The 
dashed lines are the interacting case with SLy5 force. The 
continuous lines correspond to several parameterizations of 
the tensor contribution following the T11-T16 chain. The 
system densities considered in the three panels are p — 0.08, 
0.16 and 0.24 fm~^ 
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p=0.16fm"^ ^ 
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FIG. 11: (Color online) The relative neutrino mean free path 
in neutron matter A/Afg as a function of the density for some 
parameterizations of the tensor terms. The incoming neutrino 
energy is £^=10 MeV. 



response [1, H^] , in particular in the S ~ 1 channel, are 
discussed in connection with the sum rules. 

We applied our results to the study of the neutrino 
mean free path. The advantage of the present frame- 
work is that it allows to describe nuclear (and neutron) 
matter equation of state and the neutrino mean free path 
simultaneously, hence in a self-consistent way. Obviously, 
before to achieve reliable description of neutrino trans- 
port phenomena in neutron stars, the calculations per- 
formed here must be generalized at finite temperature 
and for asymmetric nuclear matter. Nonetheless, already 
at this oversimplified level (pure neutron matter at zero 
temperature) we have shown the strong dependence of 
the neutrino mean free path on the tensor term param- 
eterizations. It represents an important reason, among 
others, to an accurate treatment of Skyrme functionals 
including tensor contribution. 
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IV. SUMMARY AND CONCLUSIONS 



In this article we have calculated the RPA response 
function for pure neutron matter considering Skyrme en- 
ergy density functionals including tensor terms. This 
article parallels with [l|, [2] where similar calculations 
were performed for the symmetric nuclear matter. As 
in previous articles divergences and instabilities of the 
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Appendix A: Particle-hole matrix elements in 
presence of a zero range tensor interaction. 

Following the notation adopted in article I and II, we 
jive in Table U the values of the particle- hole residual in- 



teraction for the tensor part of the functional in terms of 
the with x = As, F, . . . , coefficients of the functional. 
This particular notation has been already discussed in ar- 
ticle II and we refer to it for detailed explanations. 



TABLE I: Contribution of the EDF tensor part to the residual interaction in terms of the coupling constants. For the 
sake of simplicity we have introduced Kij = [{ki2)i(ki2)j], where (^12)^^ is defined in Eq.(9) of article I. The term SggiSsi is 
implicit everywhere. 





M' = 1 


M' = 


M' = -1 


M = 1 




-4B^K_i,o 


-4B^K_i,_i 


M = 


4B^Ko,i 


-2q' {B^ -4B^' +4B^' + B^) 
+4: {B'^ + B^) Ko,o 
-8B^Ki,_i 


4B^]K_i.o 


M = -1 


-4B^Ki,i 


-4B^]Ki,o 


-2g2 (B^+4B^=) 
+4B^Ko,o 
-4 (2B^ + Ki,_i 



Appendix B: Response functions 

This appendix contains the explicit expressions for the response functions for pure neutron matter. Since the isospin 
is no longer a relevant quantum number, each channel is denoted as (1, M) for 5 = 1 or only (0) for 5* = 0. 
We have 

• for the 5 = channel 



Xhf -i fJ7(o,o) , „,(o) q „, 2 

= 1 - 'xa + ( y Xo - 2kpX2 

Xr'pa ^ 



(0)12^,4 



X2 - XoXi + 



xl 



-q Xo 



2X0 



(0) 



371- 



(Bl) 



• for the S — I channels 
Xhf 

(1,0) 
Xrpa 



klm*C^ 



+ Xo 



W. 



(1) 



Xo - 2kpX2 H TT^ iXo - X2) 



37r2 



k%xl - kpXoXi + m*'^^'''Xo 



*2, ,2, 2 



jL3 

77 



2^*2^2 {wi'H2cn 



-Xo 



(B2) 
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and 



XHF 
(l.±l) 

Xrpa 



-Xa 



■2,, r^F'^*^% 



37r2 



(Xo - X2) 



hp 



Xi - XOXA + 



Xo 



q Xo 



2xc 



.(1) 
' 2 



37r2 



±i)/2) 



m * fe|. 



(1) 



X(i.±i)/2 



The coefScients W>^j 2 defined as 



(0) 



2(Co^° + Cf) + (2 + 7)(1 + 7)(C/;''^ + CDp^ - (^2^0^^" + 2Cf " + hjl 



'-.CI 1 



2 ^ 



(0) 



2 2 



(1) 



2Q 



As 



2C{ 



As 



We have also defined the w'^'^^ and X^^'^^ coefficients as 



(0,0) 



(B3) 



(B4) 



(B5) 



(1.0) 









rn J J 







1 „2 



(i,±i) 



W^^^ + 2g4 [C^J] 



(/32 - /33) 



1 + g2 (/32 - /33) 



(0) 



- 2C^ 



16 



2q2 [C^ 



2q2 [C^ 



XO - 

(/32 - /33) 



(1) 



1 + g2 (/32 - /33) [W^. 
1 + g2 (/32 - /Ja) M^i'^ 



X2 + kpXi 



, (B6) 



(B7) 

(B8) 
(B9) 



Appendix C: Landau approximation 



Since we have no isospin quantum number, the ph interaction is reduced to three terms. As done in article II we 
take the limit q and qi_2 of the second functional derivative and we obtain 



2kp [1 - cos ( 



(CI) 
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where the symbols Sab has been defined in refs.[43|, \A4 

calculated from the wj;^'' coefficients defined in Eq. (jB4[) taking the limit of W,^''^ for q — > 
So the Landau parameters can be written as 



(s) 

The various coefficients with i 

(S) 



1, 2 can be easily 



(1) 

2,L 



where Nq ^ 



g^i^^ is the normalization factor and g = 2 is the degeneracy in PNM. 
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